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ABSTRACT: The influence of environmental factors on biotic responses to nutrients was examined in three 
diverse agricultural regions of the United States. Seventy wadeable sites were selected along an agricultural 
land use gradient while minimizing natural variation within each region. Nutrients, habitat, algae, macroinver- 
tebrates, and macrophyte cover were sampled during a single summer low-flow period in 2006 or 2007. Continu- 
ous stream stage and water temperature were collected at each site for 30 days prior to sampling. Wide 
ranges of concentrations were found for total nitrogen (TN) (0.07-9.61 mg/1) and total phosphorus (TP) 
(<0. 004-0. 361 mg/1), but biotic responses including periphytic and sestonic chlorophyll a (RCHL and SCHL, 
respectively), and percent of stream bed with aquatic macrophyte (AQM) growth were not strongly related to 
concentrations of TN or TP. Pearson's coefficient of determination (R ) for nutrients and biotic measures across 
all sites ranged from 0.08 to 0.32 and generally were not higher within each region. The biotic measures (RCHL, 
SCHL, and AQM) were combined in an index to evaluate eutrophic status across sites that could have different 
biotic responses to nutrient enrichment. Stepwise multiple regression identified TN, percent canopy, median rif- 
fle depth, and daily percent change in stage as significant factors for the eutrophic index (R 2 = 0.50, p < 0.001). 
A TN threshold of 0.48 mg/1 was identified where eutrophic index scores became less responsive to increasing 
TN concentrations, for all sites. Multiple plant growth indicators should be used when evaluating eutrophica- 
tion, especially when streams contain an abundance of macrophytes. 
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INTRODUCTION to aquatic ecosystems. Agricultural lands are com- 

monly associated with high loading of non-point 
source nutrients to surface waters. The United States 
Nutrient enrichment from anthropogenic sources (U.S.) Geological Survey's (USGS) National Water 
has long been recognized as a significant impairment Quality Assessment (NAWQA) Program reported that 
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streams draining agricultural watersheds transported 
high concentrations of phosphorus and nitrogen 
(Mueller and Spahr, 2006). The U.S. Environmental 
Protection Agency's (USEPA) nationwide assessment 
on wadeable streams determined nitrogen and phos- 
phorus concentrations were elevated above reference 
(or least disturbed) sites in 53 and 47% of the 
streams studied, respectively (USEPA, 2006). 

Nutrient enrichment of streams from man-made 
sources can often cause eutrophication (i.e., increased 
primary productivity), which may result in excessive 
plant growth followed by severe diurnal variations in 
dissolved oxygen, reduction in available instream 
habitat, and a decrease in overall value for human 
uses (Welch, 1980; Dodds and Welch, 2000; Dodds 
etal., 2009). These eutrophic factors can lead to 
reductions in aquatic biodiversity that favors invasive 
species over native species (Wang et al., 2007; Maret 
et al., 2008). These concerns combined with the diffi- 
culties in effectively controlling nutrients have led 
efforts by the USEPA to develop regional nutrient cri- 
teria using an ecoregional framework to help prevent 
nuisance growth of algae in streams (USEPA, 2000a). 
Two common biotic indicators of stream nutrient 
enrichment include chlorophyll a in periphyton 
(RCHL) and sestonic (SCHL) algae. Percent coverage 
of stream bottoms by aquatic macrophytes (AQM) 
including aquatic angiosperms, bryophytes, and fila- 
mentous benthic algae (Allan, 1995) has received 
much less attention (Hering et al. , 2006), even though 
excessive macrophyte growth can have major effects 
on beneficial uses of streams (Carpenter and Lodge, 
1986; Suplee et al., 2009). 

Numerous studies have evaluated the influence of 
nutrients in experimental and natural stream chan- 
nels and the concentrations or thresholds that stimu- 
late plant growth (reviewed by Borchardt, 1996). 
Empirical regression models that predict algal bio- 
mass as a function of nutrient concentrations are often 
used to establish nutrient concentrations to protect 
streams (Lohman etal., 1992; Dodds etal., 2002; 
Stevenson etal., 2006). However, the way nutrients 
interact with physical controls in regulating algal 
biomass in streams is not well understood. Nutrient 
concentrations often explain only a small to moderate 
amount of the variation in algal biomass (Munn et al. , 
2010) and complex interactions of environmental fac- 
tors including turbidity, temperature, light, substrate, 
grazer intensity, and hydrologic regime can impact 
accrual or loss of algal biomass (Lamberti and Resh, 
1983; Munn et al., 1989; Biggs, 1995, 1996a; Figueroa- 
Nieves etal., 2006). Munn etal. (2010) determined 
the inclusion of the habitat variables temperature, 
canopy, velocity, slope, and base-flow index with nutri- 
ent concentrations improved the prediction of algal 
biomass in agricultural streams. An uneven spatial 



distribution of algae on a stream bed and variation in 
growth and biomass accrual over time may contribute 
to poor correspondence between algal biomass and 
nutrients in field studies (Morin and Cattaneo, 1992; 
Stevenson et al., 2006), but biomass and nutrient con- 
centrations can also vary diurnally and with daily or 
seasonal weather-related hydrologic events (Perkins 
and Jones, 1994; Figueroa-Nieves et al., 2006). 

Understanding how algal biomass responds to 
streamflow, habitat, and nutrient enrichment is criti- 
cal in the development of effective management strate- 
gies. Thus, it is important to examine the direct factors 
among different landscapes responsible for influencing 
stream eutrophication. For example, Munn et al. 
(2009) found local factors (i.e., substrate and water 
temperature) were more important in controlling 
macroinvertebrates than landscape features such as 
ecoregions. Some studies that include sites with a 
gradient of nutrient concentrations have been effective 
at identifying factors responsible for eutrophication 
(Lohman et al., 1992; Munn et al., 2002; Porter et al., 
2008) and correlative or regression responses that 
model algal biomass and biotic responses to environ- 
mental factors can improve understanding of algal- 
nutrient relations in streams (Porter et al., 2008). 

Streamflow is an important variable that limits 
flora and fauna and strongly influences the functional 
status of stream ecosystems (Poff, 1997) and may 
explain in part why strong nutrient-biomass relations 
have not been found in many stream studies (Lohman 
et al., 1992). Biggs and Close (1989) recommend that 
algal biomass data always be viewed within the con- 
text of antecedent streamflow conditions. High stream- 
flow events can dislodge algae or can be beneficial to 
algal communities if grazing herbivores are reduced 
(Koetsier, 2005) whereas low streamflows can cause 
desiccation of algal communities (Stanley et al., 1997). 

Flow regimes also directly influence the availabil- 
ity of nutrients and dissolved oxygen for algal growth 
and can influence water temperature and turbidity 
that directly affects water clarity. However, continu- 
ously measured streamflow data are often unavail- 
able for many stream locations, and can be expensive 
to collect. However, McMahon et al. (2003) showed 
that for certain hydrologic aspects (e.g., measure- 
ments of duration and relative change in flow), stage 
data were as useful as streamflow based metrics. 
Thus, this study investigated if stage data could char- 
acterize hydrologic conditions prior to sampling to 
better evaluate nutrient-algal growth as it relates to 
development of nutrient criteria. 

The objectives of this study were to: (1) compare 
the responses of aquatic plant growth (i.e., RCHL bio- 
mass, SCHL concentration, and percent AQM) to 
nutrient concentrations among streams in three 
diverse agricultural regions of the U.S.; (2) evaluate 
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whether proximate environmental factors influence 
biological responses to nutrients; and (3) examine 
changes in biotic responses along nutrient gradients 
for thresholds that may lead to nuisance plant 
growth. Understanding the influence of streamflow 
and habitat on nutrient-algal relationships among 
different ecoregions would help establish the range of 
expected conditions and the factors that regulate 
those ranges. USEPA (2000a) recommends using 
information about habitat and streamflow conditions, 
but offers no specific guidelines on integrating these 
variables into nutrient criteria development. This 
paper strengthens current approaches to assess 
stream eutrophication and to develop nutrient crite- 
ria for major agricultural regions of the U.S. 



METHODS 



Study Area Description 

This study was conducted in three regions (or study 
areas) characterized by a gradient of agricultural land 
use practices ranging from minimal to intense (Fig- 
ure 1, Table 1). These included: the Ozark Highlands 
(OZRK) in Arkansas, Missouri, and Oklahoma; Upper 
Mississippi (UMIS) in Minnesota and Wisconsin; and 
Upper Snake River (USNK) in Idaho and Nevada. 
Field activities were conducted in 2006 at the OZRK, 
and in 2007 at the UMIS and USNK. The USEPA 



(2000b) has developed ecoregional nutrient criteria 
that were intended as starting points for states and 
tribes to develop more refined and locally relevant 
nutrient criteria. For the three study areas in this 
study, total nitrogen (TN) and total phosphorus (TP) 
ecoregional criteria range from 0.31 to 0.54 mg/1 and 
0.010 to 0.033 mg/1, respectively (Table 1). 

All three regions had significant agricultural lands 
ranging from an average of 26% in the USNK to 49% 
in the UMIS. Livestock farming is common among all 
three regions. Long-term average precipitation for 
the OZRK, UMIS, and USNK are about 123, 78, and 
45 cm/year, respectively. The eastern regions UMIS 
and OZRK are in biomes with increased humidity 
with agriculture relying primarily on natural rainfall 
and smaller drainage basins (142-249 km 2 ). The 
OZRK region is entirely within the Ozark Highlands 
aggregate ecoregion and the UMIS region lies within 
the Glaciated Dairy Region aggregate ecoregion 
(Table 1). The western region (USNK) is located in 
the Xeric West aggregate ecoregion, relies more on 
irrigation practices, has larger average drainage 
basins, and higher elevation. High streamflows in 
this region are associated with spring snowmelt 
rather than rainfall events. Descriptions for aggre- 
gate ecoregions are found at USEPA (2000b). 



Site Selection and Study Design 

A combination of a geographic information system 
and field reconnaissance methods were used to select 
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FIGURE 1. Location of the Streams Sampled Within the Ozark Highlands (OZRK), Upper Mississippi (UMIS ), and Upper Snake (USNK) 
Regions (gray-shaded areas) Sampled by U.S. Geological Survey, 2006-2007. See Appendix 1 for names and locations of sample sites. 
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TABLE 1. Summary of Dominant Physical Characteristics (values are 
means with ranges) and USEPA (2000b) Ecoregional Nutrient Criteria. 



Parameter 


A'7TII7' / art i 

OZK.K (ra = 22) 


UMIb (ti = 18) 


TTOTWTT7" t nn\ 

UaJNK (n = SO) 


Hjcoregion level 111 


Ozark Highlands 


JN orth central Hardwood r orest 


onake Iviver Jr lam JN orthern 








Basin and Range 


Aggregate nutrient ecoregion 


Central and Eastern 


Glaciated Dairy xtegion 


Xeric West 




Forested Uplands 






Climate 


Temperate Highlands 


Humid Plains 


Arid Intermontane 


Site elevation (m) 


258 (119-388) 


320 (213-432) 


1,315 (776-1,914) 


Basin size (km 2 ) 


142 (75-255) 


249 (31-634) 


815 (0.22-5,225) 


Agricultural land (%) 


39.4 (0.5-86.9) 


48.5 (5.1-94.8) 


25.7 (<0.1-95.9) 


Precipitation (cm/year) 


123 (118-130) 


78 (69-88) 


45 (22-75) 


Ecoregional nutrient criteria 








Total nitrogen (mg/1) 


0.31 


0.54 


0.38 


Total phosphorus (mg/1) 


0.01 


0.033 


0.022 



Notes: Data summaries from U.S. Geological Survey Geographic Information (GIS) Sources. See Brightbill and Munn (2008) for a summary 
of methods used to derive GIS summaries. Ecoregion level III names from Omernik (1987); aggregate nutrient ecoregions from USEPA 
(2000a). 



the 70 study sites. Site selection attempted to identify 
study sites within a single ecoregion having similar 
hydrologic landscape, soils, climate, land use, and 
biota (Omernik, 1987). Wadeable streams with a gra- 
dient of agricultural land use and range of nutrient 
conditions were selected. Independent basins were 
selected and the National Hydrologic Dataset used to 
verify the locations of streams (Brightbill and Munn, 
2008). The smallest basins were indicative of large 
springs located in the USNK. Reconnaissance surveys 
were completed for each candidate stream to evaluate 
access, habitat conditions, and stream size. When 
existing information was not available, field measure- 
ments of nutrients, dissolved oxygen, pH, specific con- 
ductance, and temperature were measured. Field 
forms were completed that documented observations 
for habitat quality and flow characteristics. Basins 
containing USGS continuous gages were targeted if 
they met the desired conditions. Information on 
nutrient concentrations were most important in the 
site selection process to insure sites captured a gradi- 
ent of nutrients found in each study area. This pro- 
cess resulted in the selection of 18 to 30 sites within 
each region that had complete datasets for analysis 
(Appendix SI). A more detailed description of the site 
selection process can be found in Brightbill and Munn 
(2008). 

Biological and habitat assessments targeted the 
growing season (July-August) when streams reached 
stable summer low-flow conditions. The focus of this 
assessment was at the reach and microhabitat (rif- 
fle) scale using standard sampling (Table 2). A 
reach was defined as a repetition of a geomorphic 
sequence (e.g., two riffles and two pools), or 20 
channel widths if repetitive units were not available 
and at least 150 m in length (Fitzpatrick et al. , 
1998). Many of the instream habitat variables 



selected for analysis follow the conceptual template 
suggested by Biggs (1996b) to regulate aquatic plant 
growth in streams (e.g., nutrients, light, streamflow, 
substrate, temperature, and grazing by macroinver- 
tebrates). 



Water Chemistry 

Nutrient samples were collected twice at each site 
with the first sample collected ca 30 days prior to the 
biological and the second sample collected concur- 
rently with biological samples. A paired t-test found 
there was no significant difference (p > 0.05) in nutri- 
ent concentrations between time periods, so a mean 
nutrient concentration for each site was used in the 
analysis. Mean differences for paired TN and TP con- 
centration (n = 70) were 0.74 and 0.022 mg/1, respec- 
tively. Nutrient samples were collected using a depth- 
and width-integrated sampling method or one to three 
mid-channel grab samples composited when flow 
velocities and/or water depths were too low to allow 
isokinetic sampling using field equipment (Shelton, 
1994). All samples were placed on ice and shipped 
within 24 hours of collection to the USGS National 
Water Quality Laboratory (NWQL) in Arvada, Colo- 
rado, for analysis. Instantaneous field measurements 
included conductivity and discharge. Turbidity was 
measured with a Hach meter model 2100P (Hach, 
Loveland, Colorado). Nutrient samples were analyzed 
for TN and TP. Alkaline persulfate digestion was used 
to determine TN (Patton and Kryskalla, 2003) and mi- 
crokjeldahl digestion was used for TP (Patton and 
Truitt, 2000). TN and TP indicate nutrient availabil- 
ity as well or better than dissolved parameters and 
are most commonly used in nutrient criteria develop- 
ment (Dodds and Welch, 2000). Method detection 
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TABLE 2. Summary of Chemical, Physical, and Biological Variables Calculated or Measured for All Sites and Each Study Area, 2006-2007. 



Study Area 

Variable Abbreviation Units OZRK (n = 22) UMIS (n = 18) USNK (n = 30) 



Chemical 



Total nitrogen 


TN 


mg/1 


0.07-4.71 


0.46-9.61 


0.11-3.91 


Total phosphorus 


TP 


mg/1 


<0.004-0.062 


0.033-0.361 


0.010-0.159 


Conductivity 


COND 


|.iS/cm 


268-529 


252-1,089 


82-1,560 


iiysical 












Turbidity 


TURB 


NTU 


<0.1-8.0 


<0.1-50.0 


0.1-31.5 


Canopy shading 


CAN 


% 


12.0-69.7 


10.2-83.7 


6.2-70.1 


Discharge 


Q 


m 3 /s 


0.029-0.461 


0.003-0.688 


0.033-3.416 


Gradient 


GRAD 


% 


0.081-0.546 


0.013-0.367 


0.027-3.200 


Daily change in stage 


D STAGE 


% 


0.1-8.6 


0.4-3.8 


<0.1-2.8 


Base flow index 


BFI 


% 


26.4-49.6 


45.6-67.4 


64.6-86.6 


High flow index 


HFI 


Unitless 


0.376-2.418 


0.682-1.771 


0.347-1.286 


Dimensionless shear stress 


DIMSHEAR 


Unitless 


0.000-0.030 


0.000-0.020 


0.000-0.040 


Water temperature 


TEMP 


C 


18.6-26.1 


14.2-25.8 


14.5-24.9 


Median substrate 1 


SUBSTR 


mm 


1-96 


9-384 


9-384 


Median depth 1 


DEPTH 


m 


0.067-0.213 


0.030-0.274 


0.08-0.50 


Median velocity 1 


VEL 


m/s 


0.07-0.81 


0.02-0.48 


0.03-0.59 


iological 












Periphyton chlorophyll a 1 


RCHL 


mg/m 2 


10.4-124.7 


8.3-172.3 


2.4-149.9 


Seston chlorophyll a 


SCHL 


Mg/1 


0.2-2.1 


0.5-36.4 


0.1-6.0 


Aquatic macrophytes 


AQM 


% 


0.0-50.0 


0.0-62.9 


0.0-90.7 


Invertebrate scrapers 1 


SCBIO 


g/m 2 


0.1-25.5 


<0. 1-7.4 


<0. 1-17.7 


Eutrophic index 


EI 


Unitless 


6.3-28.3 


5.9-44.7 


1.2-37.6 



1 Variable represents microhabitat (riffle) within reach. 



levels for TN and TP were 0.03 and 0.004 mg/1, 
respectively. Censored values were rare, with only 
one occurrence for TP, which was assigned a value of 
one-half the method detection level prior to analysis. 
More details about sampling or laboratory methods 
can be found in Brightbill and Munn (2008). 

Biological Samples 

A Slack sampler (500 micron mesh, 0.25 m 2 sample 
area, total area of 1.25 m 2 per sample) was used to 
collect semi-quantitative invertebrate samples from 
five riffle locations throughout the reach. Sample pro- 
cessing involved elutriation to remove inorganic deb- 
ris, compositing subsamples, and preservation with 
10% formalin (Moulton etal, 2002). Macroinverte- 
brates were sorted using standard 500 count proce- 
dures, identified to the lowest practical taxonomic 
level, and biomass (g/m 2 ) determined (Moulton et al., 
2000). Biomass estimated were determined by drying 
at 60°C to a constant weight. Large-rare taxa (e.g., 
crayfish) were not included in the biomass determina- 
tions. The biomass of scrapers was determined based 
on the ratio of macroinvertebrates that were scrapers. 
Trophic groupings were assigned using the USGS 
Invertebrate Data Analysis System (Cuffney, 2003). 

Chlorophyll a samples were collected from coarse 
substrate (gravel or larger) and the water column. 



For RCHL samples, five rocks were collected in close 
proximity to each of the invertebrate samples and the 
algae were scraped from a consistent area on each 
rock and 25 subsamples composited. Water column 
samples for SCHL were collected by methods used for 
nutrient sampling. RCHL and SCHL samples were 
filtered onto a Whatman 47-mm glass fiber filter and 
analyzed for chlorophyll a and pheophytin biomass 
(Moulton et al. , 2002). Algal biomass samples were 
shipped on dry ice to the NWQL for analysis. Chloro- 
phyll a and pheophytin were determined using an 
adaptation of EPA method 445.0 (Arar and Colling, 
1997), with a detection limit of 0. 1 u.g/1. 

Visual estimates of percent AQM (including macro- 
algae with filaments >2 cm) coverage were made at 
five points across each of 11 equidistant transects 
along the study reach. A mean percentage for the 55 
observations was calculated to represent the reach. 
Forty percent coverage of the stream bottom with 
AQM was selected to characterize eutrophic condi- 
tions with excessive plant growth. This was derived 
by taking the midpoint of 20-60% coverage reported 
by Suplee et al. (2009) as undesirable growth in Mon- 
tana streams. Similarly, Chambers et al. (1999) 
reported levels of 10 to 50% may be considered nui- 
sance growth in Canadian streams. We also identified 
a group of sites that appeared to have excessive 
AQM growth covering more than 40% of the stream 
bottom. 
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Habitat 

Physical habitat was assessed at the stream reach 
and microhabitat (riffle) scale. A total of 11 equidis- 
tant transects oriented perpendicular to streamflow 
were established throughout the reach, with channel 
width (m) measured at each transect. Mean percent 
canopy cover was measured using a densiometer with 
measurements made at the center channel and each 
bank edge along each transect. Measurements were 
averaged for each transect and then combined for an 
overall mean percent canopy cover for each reach. 
Reach gradient was determined using a surveyor's 
level with gradient calculated at the water-surface. 
Characterization of riffle microhabitats consisted of 
measurements of near-bed velocity (ca 2 cm from 
stream bottom), depth, and substrate size at each 
location where algae and macroinvertebrates were 
sampled. A modified Wolman pebble count was com- 
pleted at riffle habitats sampled (Wolman, 1954). This 
consisted of selecting and measuring substrate parti- 
cles at each of five locations where macroinverte- 
brates and algae were sampled using a square PVC 
frame with 20 equally spaced locations around the 
frame. A total of 100 particles were measured to 
characterize the substrate size into 10 substrate size 
classes for the riffle habitats sampled. Median values 
of velocity, depth, and substrate were calculated to 
represent microhabitat measurements. Additional 
detail on methods used to collect habitat data can be 
found in Fitzpatrick et al. (1998). 

Continuous Stage and Water Temperature 

Continuous stream-stage data and water tempera- 
ture were collected at existing USGS gages (23 sites) 
or vented submersible pressure transducers (47 sites) 
with a 15-min to 1-hour time step. Stage and tempera- 
ture data from sites within each region were pro- 
cessed to 1-hour intervals and used to derive daily 
statistics using a common period of record of 30 days 
prior to the biological sampling (McMahon et al. , 
2003). This period of record has been suggested to be 
adequate for biomass accrual in streams (Cattaneo 
and Amireault, 1992; Lohman et al., 1992; Snyder 
et al., 2002) and it allows for maximum biomass 
accrual, but avoids subsequent sloughing. 

Stage data, which have arbitrary datum, were 
adjusted using the mean bed elevation so that the daily 
values of stage provided a depth time series. Mean 
daily percent change in stage for the 30 days prior to 
the ecological assessment was calculated and normal- 
ized by median stage for this same 30-day period. 

A surveyed channel cross-section at the location of 
the transducer was used to calculate cross-sectional 



area. Discharge was measured at the time of sam- 
pling at each site using methods described in Rantz 
et al. (1982) or using established rating curves at co- 
located USGS gages. Daily-stage time-series data for 
each site were inspected for continuity and extreme 
high- or low-flow events. Dimensionless shear stress 
(DIMSHEAR), which provides measure of the force of 
streamflow relative to the size of bed material, was 
calculated for each site to evaluate bed stability: 

DIMSHEAR= MAXSTAGE x FRICTION SLOPE 

1.7xSUBSTR ' 



where MAXSTAGE was the maximum daily stage for 
the period of record and FRICTION SLOPE was esti- 
mated from the vertical velocity gradient and bed sub- 
strate (SUBSTR). A threshold of DIMSHEAR > 0.045 
was used to assess whether high flows were likely to 
have entrained bed material during the 30 days prior 
to sampling (Konrad et al., 2002). 

Streamflow records from long-term USGS gages 
were used to estimate annual maximum daily 
streamflow at each site for each of the 10 years prior 
to biological sampling (Dave Wolock, USGS, January 
15, 2009, personal communication). A high-flow 
index, representing the magnitude of high flow dur- 
ing the year of sampling relative to median annual 
high flow, was calculated as maximum daily stream- 
flow for the year of sampling at each site normalized 
by median annual maximum streamflow for the 10- 
year period. Base-flow index, the component of 
streamflow that can be attributed to groundwater dis- 
charge into streams, was estimated for each site from 
the national base-flow index 1-km resolution dataset 
developed by Wolock (2003). 

Data Analyses 

All statistics were calculated using either Systat® 
version 11.0 (Wilkinson, 2004) or R (R Development 
Core Team, 2005). Scatter plots were examined for 
outliers and spurious correlations. Correlation coeffi- 
cients (Spearman's rho) among variables were exam- 
ined to identify and reduce redundancy. Variables 
with skewed distributions were logio or square-root 
transformed to normalize their distributions. Follow- 
ing transformation of the data, ANOVA was used to 
determine if there were significant differences in 
variables among regions, followed by Tukey multiple 
comparisons 

Bivariate linear regression models were developed 
to examine the influence of nutrients on eutrophic 
indicators (RCHL, SCHL, and AQM). Piecewise 
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regression or segmented regression (Toms and Les- 
perance, 2003) was used to identify breakpoints or 
thresholds in biotic responses and nutrient concentra- 
tions. This approach to modeling data identifies 
regression changes at one or more points along the 
range of the independent variable. A LOWESS 
(LOcally WEighted Scatterplot Smoothing) technique, 
a robust nonparametric description of data patterns 
(Helsel and Hirsch, 2002), was initially used to exam- 
ine the number and location of the breakpoints. The 
influence of habitat and biotic factors on nutrient- 
eutrophication response was assessed using multiple 
regression with the RCHL, SCHL, or AQM as the 
dependent variable and TN, TP, and selected habitat 
variables as the potential independent variables (see 
Table 2). Conductivity was evaluated as an explana- 
tory variable as it has been shown to be a surrogate 
for nutrients (Biggs, 1995; Munn et al., 2002). 

An exhaustive stepwise search identified the mul- 
tivariate models with the lowest Akaike's informa- 
tion criterion (Helsel and Hirsch, 2002) where all 
independent variables had statistically significant 
coefficients (p < 0.05 that the coefficient = 0) and 
multiple colinearity was not an issue (variance infla- 
tion factor <2 for all explanatory variables). Three 
binary variables indicating regional membership 
were created to determine if region-specific models 
were warranted. 

A eutrophic index (EI) was calculated using three 
metrics that reflect instream primary productivity 
(RCHL, SCHL, and AQM) to nutrients expressing 
eutrophication. Each metric was scaled to an 
expected low and high for all sites so the range of 
scoring approximated 0-100. This is similar to the 
scoring system outlined in Hughes et al. (1998): 



El 



RCHL/2 + SCHL x 2 + AQM 
3 



(2) 



This empirical EI is not intended to be transfera- 
ble to other investigations because it is based on the 
observed range of biotic responses at our sites. Multi- 
ple linear regression, as described above, was used to 
develop a model expressing the influence of nutrients, 
habitat, and biotic factors. 



RESULTS 



Physical and Chemical Characteristics 

Daily-stage records from 30 days prior to biological 
sampling indicated no high flows capable of substan- 
tial bed material entrainment (DIMSHEAR < 0.045) 



or extreme low flows in reach habitat sampled. A few 
sites had high values (DIMSHEAR = 0.040), but were 
stable at this flow. At these sites, the reaches had 
boulder steps with very high shear stress. As a 
result, the reach average shear stress calculations 
were higher than and unrepresentative of shear 
stress at the location of bed material and biological 
sampling. Daily change in stage was generally low at 
the sites (median of about 1%), though it did range 
up to about 9% (Table 2). Mean high-flow index val- 
ues for each region were near 1.0 indicating the year 
of sampling was generally similar to the long-term 
10-year period prior to sampling. Mean base-flow 
index values among study regions were significantly 
different, with the USNK sites having the largest 
groundwater contributions to total streamflow (76%). 
There were no statistical differences among regions 
for conductivity, turbidity, and canopy shading. 

The sites had a gradient of TN and TP concentra- 
tions within and among regions (Figures 2a and 2b) 
and represented a wide trophic status from low- to 
high-nutrient conditions. The gradient and range for 
TP concentrations was the weakest among OZRK 
sites. Concentrations of TN ranged from 0.07 to 
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FIGURE 2. Cumulative Distributions of (a) Total Nitrogen (TN) 
and (b) Total Phosphorus (TP) Across All Study Sites, n = 70. 
Dashed lines indicate nutrient (TN and TP) trophic boundaries 
suggested by Dodds et al. (1998). L, low enrichment; M, moderate 
enrichment; H, high enrichment. 
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9.61 mg/1 for all sites (Table 2). Concentrations of TP 
ranged from <0.004 to 0.361 mg/1, with the UMIS 
having significantly greater mean concentration 
(0.128 mg/1) and the OZRK significantly lower mean 
concentrations (0.024 mg/1) than other regions. There 
was a positive relation between TN and TP (r = 0.55, 
p < 0.05) (Table 3). Conductivity was also positively 
related to TN (r = 0.43, p < 0.05). 

The classification scheme of Dodds et al. (1998) 
indicated that about 26-30% of all sites were highly 
enriched based on the TP and TN concentrations, 
respectively (Figure 3). The UMIS region had the 
highest percentage of enriched sites with 50 and 
56% of the sites exceeding TP and TN concentra- 
tions, respectively. The OZRK region had only about 
23% enriched sites based on TN concentrations and 
no sites exceeding the highly enriched boundary for 
TP. 



Biological Measures of Eutrophication 

Biomass of RCHL ranged from 2.4 to 172.3 mg/m 2 
for all sites with no significant difference between 
regions (Table 2). The UMIS had the lowest mean 
RCHL biomass (46.6 mg/m 2 ) among regions despite 
the highest mean TN and TP concentrations among 
regions. Concentration of SCHL ranged from 0.1 to 
36.4 ng/1, with the UMIS mean (7.6 ng/1) significantly 
higher than the other regions. The AQM were highly 
variable and ranged from 0 to 91%, with no signifi- 
cant differences in regions. The highest visual per- 
centages for AQM were found in the USNK region 
where Potomogeton sp. or Cladophora sp. often com- 
prised the majority of the streambed coverage. A posi- 
tive relation was found between AQM and median 
depth (r = 0.55, p < 0.05) (Table 3). 

High values of AQM and RCHL were inversely 
related with sites forming a wedge-shaped scatterplot 
(Figure 4). There was a group of sites with AQM that 
approached or exceeded 40% bottom coverage within 
the study reach that indicated excessive plant 
growth. However, many of these same sites had rela- 
tively low (<50 mg/m 2 ) RCHL biomass. 

The invertebrate scraper biomass ranged from 
<0.1 to 26 g/m 2 for all sites (Table 2). Molluscs 
exceeded 3 g/m 2 at 16, 4, and 8 sites in the OZRK, 
UMIS, and USNK study units, respectively. The 
greatest mean values (5.7 g/m 2 ) were found for the 
OZRK sites, followed by USNK (2.5 g/m 2 ) and UMIS 
(1.0 g/m ). The largest contributor to the biomass in 
the OZRK was Elimia potosiensis, a native snail. 
The USNK sites with the greatest biomass were a 
result of Potamopyrgus antipodarum, an invasive 
snail. There were no statistically significant correla- 
tions between invertebrate scraper biomass and 



O 



II 

% o 
o 

O © 

a t> 
.2 ii 

-+J 

h 2 

u ni 

o <u 

a 

« 2 

-a M 

CU 

p, Z. 
CO -a 
d 

CO ® 

H CO 

j a 



o 



— 

o 



o 



rH IS 



O 
I 



rH OS LQ LO 
rH O rH O 

dodo 



CO LO CO O 

rH O rH CN rH 

<6 d <6 d <6 

I I 



o o o o o 
I I I 



o o o o o o o 
I I 



























rH 


CM 


CO 


o 


CO 




o 




























o 


o 


o 


rH 


rH 




CN 


CO 


w 
























d 


d 


o 


d 


o 


O 


d 


o 
































I 
































rH 


rH 


t> 


cc 


as 






LO 


CM 
























o 


CO 


CM 


CO 


o 


• 




o 


o 


H 






















d 


d 


d 




d 






d 




H 




























1 






1 


] 




— 


























CO 


as 


CO 


rH 


CM 


OS 


as 






















CN 


rH 


rH 


CM 


CN 


o 


rH 


rH 


o 


o 






















d 


d 


d 


d 


d 


d 


O 


d 


d 


d 




































I 


] 


I 


r— 1 








































Q 


























































CN 


o 


LO 


rH 


LO 


o 


CM 


CN 


CO 


o 


X 


Of 


















o 
d 


-0.4 


-0.0 


-0.2 


-0.0 


-# 
d 


-0.0 


rH 

d 


CO 

d 


-0.1 


rH 

d 


















CN 


LO 


o 


CM 


LO 


LO 


CD 


rH 


CO 


rH 


CO 


CM 


— 
















o 


I- 


CN 


o 


o 


CM 


CO 


rH 


O 


CM 


rH 


CO 


















d 


d 


d 


d 


d 


d 


d 


d 


d 


d 
I 


d 


d 
I 


s 














LO 


CM 


CN 


00 


CM 


CO 


CD 


CN 




CN 


10 


CD 


as 
















rH 




rH 


CN 


o 


o 


o 


CO 


o 




LO 


rH 


CM 
















d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


— 














I 




1 




I 


I 


I 






























rH 


CO 


CD 


CO 


o 


as 


CN 


o 


rH 


X 


CO 


CD 














o 


o 


o 


o 


O 


rH 


o 


o 


CO 


CM 


CN 




LO 


o 


5 












d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


H 












I 


I 






! 






I 




I 




I 
















03 


CN 


rH 


CD 


LO 


rH 


CM 




CM 


as 


as 


o 


LO 




o 












o 


O 


CM 


o 


o 


CN 


rH 


rH 


rH 


o 


o 


CN 


rH 


rH 


CO 


O 










d 
I 


d 


d 


d 

I 


d 


d 

1 


d 
I 


d 


d 


d 


d 


d 

! 


d 
1 


d 


d 
I 


« 








































H 








rH 


oo 


CN 


CD 


rH 


o 


t> 


as 


rH 


t> 


t> 


LO 


CN 


CM 


CM 


t> 


CB 








o 


CN 


o 


rH 


o 


o 


rH 


rH 


CO 


o 


rH 


o 


O 


O 




rH 


CQ 








d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


P 










I 










1 


1 


1 






I 


I 






1 


m 














































t> 


o 




CO 


as 


CN 


rH 


rH 


o 


rH 




CO 


[> 


o 




CO 


CM 








10 


CN 


rH 


o 


o 


o 


rH 


CO 


rH 


rH 


o 


CN 


o 


rH 


O 


CO 


CM 








d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 












1 










] 


I 


! 










I 


1 


I 






CO 


rH 


CO 




rH 


rH 


rH 


0) 


CD 


CO 




rH 


t> 


CM 


CN 


CK 


O 


CO 


Ph 




CN 




rH 


10 


o 


CN 


o 


o 


rH 


CM 


rH 


CM 


CN 


o 


O 


O 


LO 


CO 


H 




d 

I 


d 


d 


d 


d 


d 
I 


d 

1 


d 
I 


d 

I 


d 


d 


d 
I 


d 


d 
I 


d 


d 


d 


d 




10 


CN 


i — i 






o 


o 


LO 


t> 


rH 


CO 


rH 


t> 


o 


o 


CN 


CD 


CO 


CO 


Z 


10 


CN 


CN 


o 


rH 


rH 


o 


o 


o 


rH 




o 


o 


o 


rH 


CN 


CN 


00 


LO 


H 


d 


d 
i 


d 


d 


d 


d 


d 


d 


d 


d 

1 


d 


d 

i 


d 
i 


d 


d 


d 


d 


d 


d 



Ph 

O co O 



Ph Ph J SgZ^H„ff|K|ffi 

HO>G^cocHOcoffipqcopH<;coH 



Journal of the American Water Resources Association 



505 



JAWRA 



Maret, Konrad, and Tranmer 



100 



Q. 

o 

1- 
3 

HI 

' 

c 

o 
o 

i_ 

o 

Q. 




FIGURE 3. The Relative Percentage of All Sites and Individual 
Regions Classified as Eutrophic (highly enriched) on the Basis of 
Total Nitrogen (TN > 1.5 mg/1), Total Phosphorus (TP > 0.75), 
Periphytic Chlorophyll a (RCHL > 70 mg/m 2 ), Seston Chlorophyll 
a (SCHL > 30 ug/1) (Dodds et al, 1998), and Aquatic Macrophyte 
(AQM) Coverage of Stream Bottom >40% (Chambers etal, 1999; 
Suplee et al, 2009). 
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FIGURE 4. Relation of Periphytic Chlorophyll a (RCHL) 
to Macrophyte (AQM) Percent Coverage of Stream 
Bottom for All Study Sites (n = 70). 



other environmental variables measured at study 
sites (Table 3). 

For all sites combined, the RCHL measurements 
resulted in the most sites being classified as eutro- 
phic (26%), followed by AQM (21%) and SCHL (3%, 
Figure 3). The USNK region had the highest number 
of sites (30%) that exceeded the AQM eutrophic 
boundary, followed by UMIS (22%) and OZRK (9%). 
The USNK also had the highest number of sites 
(37%) classified as eutrophic based on RCHL mea- 



surements, followed by OZRK (23%) and UMIS 
(11%). With the exception of the UMIS, other regions 
had no sites classified as eutrophic based on SCHL 
measurements . 



Nutrient-Biotic Response Models 

Variability in the three biotic responses to nutrient 
concentrations was high in all regions (Figure 5). 
Bivariate regression analyses on all sites and individ- 
ual regions produced 9 out of 24 significant models 
among TN and TP and the three biotic response vari- 
ables (Table 4). The amount of variation explained by 
these significant models was 8-32% (Figures 5a and 
5d). TN was significant (R 2 = 0.08, p < 0.018) in pre- 
dicting RCHL for all sites (Figure 5a, Table 4). For 
all sites combined, TN and TP were both significant 
in predicting SCHL (R 2 = 0.21 and 0.32, p < 0.001) 
(Figure 5c and 5d). Because many studies have found 
significant positive correlations of conductivity with 
nutrient concentrations (Biggs, 1995; Carpenter and 
Waite, 2000), TN and TP were substituted with 
conductivity. However, this did not improve model 
prediction for the three biotic variables (R 2 = 0.02- 
0.09). 

Generally, bivariate models for specific regions did 
not explain significantly more variation in biotic 
response variables than for all sites combined. There 
were no significant models identified for the OZRK 
(Table 4). The lack of a strong TP gradient among 
OZRK sites could be a partial explanation for this. 
Significant models were identified for TN and TP in 
predicting SCHL (R 2 = 0.39, p = 0.005 and R 2 = 0.32, 
p = 0.015) for UMIS. TP was also identified as a sig- 
nificant predictor of AQM in the UMIS (R 2 = 0.26, 
p = 0.030). However, this was an anomalous relation 
compared to other significant models, where TP con- 
centration was inversely related to AQM (Figure 5f, 
Table 4). In the USNK, TN was significant in pre- 
dicting AQM (R 2 = 0.34, p < 0.001) and RCHL 
(R 2 = 0.14, p < 0.043). A significant model also iden- 
tified TP as important in predicting SCHL 
(R 2 = 0.14, p < 0.043). 

Munn et al. (2010) found an improvement in 
nutrients and algal biomass relations for open can- 
opy sites vs. shaded sites. However, bivariate 
regression model results from this study used only 
sites with <50% canopy (n = 48) and generally did 
not improve model performance (data not shown). 
There were no significant models found for RCHL 
or AQM and nutrients. However, models for SCHL 
and nutrients improved slightly for TN (R 2 = 0.33, 
p < 0.001) and TP (R 2 = 0.41, p < 0.001) (see Table 4 
for comparisons). In addition, there were no signifi- 
cant differences between means in RCHL, SCHL, 
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FIGURE 5. Bivariate Plots of the Biotic Response Variables Periphytic (a and b) Chlorophyll a (RCHL), (c and d) Sestonic 
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and USNK (dotted line). Regression equations of all sites combined and individual regions can be found in Table 4. 



and AQM for open (<50%) vs. closed (>50%) canopy 
sites. 

Models for the Eutrophic Index 

The majority of variance in the EI was explained 
by either TN (R 2 = 0.29, p < 0.001) or TP (R 2 = 0.14, 



p = 0.002) (Figures 6a and 6b), although a significant 
piecewise regression (p = 0.039) was found between 
TN and EI scores (Figure 7) with a breakpoint 
(threshold) at 0.48 mg/1 (95% confidence interval of 
0.22-0.75 mg/1). This threshold was the point where 
EI scores became less responsive to increasing TN 
concentrations. A nonparametric LOWESS smooth 
also identified a similar TN threshold. A piecewise 
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TABLE 4. Bivariate Regression Models for Concentrations of 
Periphyton Biomass (RCHL), Seston Biomass (SCHL), and Aquatic 
Macrophytes and/or Macroalgae (AQM) as a Function of Total 
Nitrogen (TN) and Total Phosphorus (TP) for All Sites Combined 
(re = 70) and for the OZRK, UMIS, and USNK Study Areas, 2006- 
2007. 



Dependent 


Independent 






Model 


Variable 


Variable 


Intercept 


K 


P 


All sites (re = 


70) 








log(RCHL) 


0.214 log(TN) 




u.uo 


U.UXo 


log(RCHL) 


-0.002 log(TP) 




a aa 
U.UU 


A QC7 


log(SCHL) 


0.474 log(TN) 


0 119 


0 91 


si\ nni 

sU.UUl 


log(SCHL) 


0.628 log(TP) 


o qq£ 


a 99 

U.oZ 


vU.UUl 


log(AQM) 


0.272 log(TN) 


o qro 


a c\a 
u.u^ 


A A8Q 
U.Uoo 


log(AQM) 


0.179 log(TP) 


1 1 Q7 


n 09 


n 9ss 


OZRK study area (re = 22) 








log(RCHL) 


0.788 log(TN) 


— l.DO / 


U. ID 


A Afi/L 


log(RCHL) 


0.056 log(TP) 


1 7CQ 


a ai 

U.Ul 


A CCO 


log(SCHL) 


0.134 log(TN) 


O 90A 


u.u^ 


U.OOZ 


log(SCHL) 


0.229 log(TP) 


o i fip; 


n i n 

U. 1U 


u. i^ty 


log(AQM) 


0.131 log(TN) 


o 809 


a ni 

U.Ul 
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log(AQM) 


0.589 log(TP) 


1 890 


n iq 

U. lo 


a ap;9 
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UMIS study 


area (re = 18) 








log(RCHL) 


0.294 log(TN) 


1.489 


0.17 


0.089 


log(RCHL) 


0.364 log(TP) 


1.930 


0.14 


0.131 


log(SCHL) 


0.900 log(TN) 


0.191 


0.39 


0.005 


log(SCHL) 


1.111 log(TP) 


1.535 


0.32 


0.015 


log(AQM) 


-0.592 log(TN) 


1.194 


0.13 


0.135 


log(AQM) 


-1.144 log(TP) 


-0.101 


0.26 


0.030 


USNK study 


area (re = 30) 








log(RCHL) 


0.421 log(TN) 


1.650 


0.14 


0.043 


log(RCHL) 


0.054 log(TP) 


1.631 


0.00 


0.845 


log(SCHL) 


0.156 log(TN) 


0.076 


0.02 


0.423 


log(SCHL) 


0.493 log(TP) 


0.719 


0.14 


0.043 


log(AQM) 


0.883 log(TN) 


1.221 


0.34 


<0.001 


log(AQM) 
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FIGURE 6. Bivariate Plots of Eutrophic Index Scores and (a) Total 
Nitrogen (TN) and (b) Total Phosphorus (TP) for All Study Sites. 
The lines indicate the best fit linear regression for all study sites. 
Dashed lines are the 95% confidence interval for the regression 
line. 



regression and breakpoint for TP and EI scores was 
not significant (p > 0.05). 

The final multivariate model for EI was a function 
of TN, canopy shading, median depth, and daily 
change in stage (adjusted R 2 = 0.50, p < 0.001) 
(Table 5) with a residual standard error of 6.6 com- 
pared to a median observed score of 16. The EI scores 
ranged from 1.2 to 44.7 for all study sites (Table 2). 
Analysis of covariance using membership in each 
region as a binary variable did not justify a separate 
regression analysis for each region. Median depth 
appeared to be a surrogate for stream size as it was 
significantly correlated with discharge (r = 0.72, 
p < 0.05) (Table 3). 

Multivariate analysis of the individual components 
of the EI identified TN and median depth as impor- 
tant predictors for RCL and AQM, respectively 
(Table 5). Several nutrient and physical parameters 
were important in predicting SCHL concentrations, 
specifically increasing with TP and conductivity and 
decreasing with median substrate size, percent 
canopy, median depth, and turbidity. 



R 2 = 0.27(n=70) 
p = 0.039 

y = 27.1 (TN) + 5.14 (TN < 0.48) 
y = 1.45 (TN) + 18.3 (TN > 0.48) 
A 




FIGURE 7. Eutrophic Index Scores as a Function of Total Nitrogen 
(TN) Concentrations for All Study Sites (re = 70). Piecewise regres- 
sion line with identified breakpoint threshold for TN of 0.48 mg/1. 
Dashed lines are the 95% confidence interval for the regression 
line. The gray line represents a nonparemetric LOWESS (Locally 
Weighted Scatterplot Smoother) fit of the data (Helsel and Hirsch, 
2002). 
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TABLE 5. Multiple Regression Models for Predicting Concentrations of Periphyton Biomass (RCHL), 
Seston Biomass (SCHL), Aquatic Macrophytes and/or Macroalgae (AQM), and Eutrophic Index (EI) 
as a Function of Chemical, Physical, and Biological Variables for All Sites, n = 70. 



Dependent 



Variable 


Predictor Variable and Coefficient 


Intercept 


Adjusted R 2 


p-Value 


RCHL 


19.04 log(TN) 


56.26 


0.05 


0.030 


SCHL 


6.55 log(TP) + -2.93 log(SUBSTR) + -0.08(CAN) 


8.47 


0.47 


<0.001 




+ -23.86(DEPTH) + 6.22 log(COND) + -3.29 log(TURB) 








AQM 


122.03(DEPTH) 


-1.05 


0.28 


<0.001 


EI 


10.41 log(TN) + -0.18(CAN) + 28.02(DEPTH) + 151.20(DSTAGE) 


18.75 


0.50 


<0.001 



Notes: Predictor variables listed were significant (p < 0.05), with most important variables listed first. See Table 2 for variable definitions. 



DISCUSSION 



Physical and Chemical Characteristics 

Evaluating nutrient-biotic responses in agricul- 
tural streams can be challenging because agricultural 
streams are known to contain some of the highest 
concentrations of nutrients from fertilizer and man- 
ure applications and commonly have altered riparian 
and instream habitats that can influence nutrient 
uptake and algal growth (Watzin and Mcintosh, 
1999; Mueller and Spahr, 2006). This can complicate 
nutrient gradient studies because of difficulty with 
locating sites with low nutrient concentrations. This 
study captured a strong TN gradient (0.07 to 
9.6 mg/1) with 43% of sites having average concentra- 
tions below concentrations (0.52 mg/1) associated with 
excessive algal growth (Dodds etal., 2006). Evaluat- 
ing these low level nutrient conditions is important to 
determine if assessment methods are sensitive 
enough to detect biotic responses in areas that are 
relatively undisturbed. 

The TN threshold of 0.48 mg/1 using the EI scores 
is in close agreement with thresholds suggested by 
Dodds et al. (2006) and Stevenson et al. (2006) (0.52 
and 0.40 mg/1, respectively) and is within the range 
of USEPA ecoregional TN criteria (0.31 to 0.54 mg/1) 
for all regions. Our identified threshold also agrees 
closely with the 75th percentile of TN concentrations 
(0.50 mg/1) measured in undeveloped watersheds 
across the U.S. (Clark etal, 2000). A TP threshold 
was not identified, even though about 36% of the 
sites sampled had average TP concentrations below 
ca 0.03 mg/1 reported to promote undesirable levels 
of algal growth (Dodds et al., 2006). 

Conductivity has been shown to be an indicator of 
water quality conditions, especially nutrient enrich- 
ment (Munn et al., 2002). Many studies have reported 
significant positive correlations of conductivity with 
nutrient concentrations, particularly in relation to 
agriculture (Biggs, 1995; Carpenter and Waite, 2000). 



We also found a significant positive relation between 
conductivity and nutrients. However, bivariate model 
performance did not improve when nutrients were 
substituted with conductivity. 

Water temperature was not important in our mod- 
els, even though other studies have found a relation 
with plant growth (Munn etal., 1989; Francoeur 
etal, 1999). The USNK sites had the lowest mean 
water temperature, primarily due to higher overall 
elevation, but plant growth was unaffected, as there 
was no significant difference among regions in RCHL 
and AQM. 

Biological Measures of Eutrophication 

We employed EI as a comprehensive measure of 
eutrophication using the three different forms of bio- 
tic responses to nutrients. To our knowledge, this 
approach has not been attempted. Our assessment 
demonstrates the need to evaluate eutrophication 
using more than one biotic response measure. Collec- 
tively, one or more eutrophic boundaries for RCHL 
(>70 mg/1), SCHL (>30 u.g/1), and AQM (>40%) were 
exceeded at 46% of sites sampled, whereas excee- 
dances based on any individual biotic measures were 
only 3-26%. The RCHL biomass sampled in this study 
represented a gradient ranging from 2.4 to 
174 mg/m 2 with about 26% of these sites considered 
highly enriched based on Dodds et al (1998). There is 
no one standard RCHL biomass considered excessive 
or nuisance level; Dodds et al's (1998) review of the 
literature found excessive periphyton chlorophyll a 
ranged from 50 to 200 mg/1. Using their recommen- 
dation of 150 mg/1, there would be only 3 (4%) of our 
sample sites exceeding this amount, whereas our 
visual assessment of percent AQM coverage 
approaching or exceeding 40% of the reach would 
indicate that 15 (21%) sample sites would have exces- 
sive plant growth (Suplee et al. , 2009). The algal 
response variables RCHL and SCHL have commonly 
been used to evaluate nutrient enrichment, whereas 
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AQM has received much less attention in monitoring 
programs. However, in many cases, RCHL or SCHL 
samples would not constitute eutrophication for these 
same sites. Indeed, high levels of both RCHL and 
AQM were not found in any streams (Figure 4). 

The nutrient-rich agricultural streams we exam- 
ined had relatively low SCHL concentrations 
(<10 p.g/1) at 93% of all sites. Only four sites in the 
UMIS had higher concentrations approaching or 
exceeding the 30 Lig/1 concentration considered to be 
highly enriched according to Dodds et al. (1998). 
Figueroa-Nieves et al. (2006) found similar low 
concentrations (2-20 p.g/1) in agricultural streams of 
Illinois. These findings suggest that measures of 
SCHL may not always be an appropriate measure 
of eutrophication in small agricultural streams. 
According to Allan (1995), the primary source of 
seston in fast-flowing streams is the sloughing of 
attached algae. However, in sluggish lowland streams 
or rivers where water mass has a longer residence 
time, conditions may be suitable for true plankton to 
colonize and reproduce. For example, in the tributar- 
ies for the Snake River sampled in USNK, the maxi- 
mum sestonic chlorophyll a concentrations was 
6 ng/1, whereas, in the free-flowing sections of the 
main-stem Snake River, concentrations approaching 
100 (ig/1 have been reported (Myers et al., 2003). 

Macroinvertebrates have been shown to control 
RCHL, especially when hydrologic disturbance is 
lacking (Jacoby, 1987; Stevenson et al, 2006). They 
found invertebrate biomass of 1-3 g/m 2 to be suffi- 
cient to reduce periphytic algae. In addition, Sumner 
and Mclntire (1982) found that in a laboratory 
stream, snails can reduce periphyton by as much as 
30%. The relatively high biomass of invertebrate 
scrapers found at many of our study sites, many of 
which contained large numbers of snails, would sug- 
gest a lack of hydrologic disturbance and potential 
herbivory control on RCHL. However, invertebrate 
scrapers were not a significant explanatory variable 
for RCHL at our sites. 



Nutrient-Habitat-Biotic Response Models 

Evaluating antecedent conditions prior to sampling 
validated using an index period approach that mini- 
mizes streamflow disturbance on plant growth. 
Continuous-stage data collected 30 days prior to 
sampling confirmed that high disturbances during 
summer index period were rare and DIMSHEAR did 
not exceed 0.045, which represents a threshold for 
bed material entrainment. High-flow disturbance was 
not a factor influencing biotic measures of stream 
eutrophication, and allowed analysis of other environ- 
mental factors on algal growth. 



Numerous lotic studies have reported that nutri- 
ents are important for algal growth, but the devel- 
opment of nutrient regression models has not been 
particularly successful. Biotic measures of stream 
eutrophication including RCHL, SCHL, and AQM 
growth were not strongly related to concentrations 
of TN or TP; R 2 values ranged from 0.02 to 0.32 
for all sites combined. These results contrast with 
Chetelat et al. (1999) who reported that nutrients 
can become a dominant variable during stable, low- 
flow conditions. Our predictive models for SCHL 
had the highest overall R z values, which may be 
partially due to autocorrelation where sestonic algal 
samples contain both chlorophyll a and nutrients 
(van Nieuwenhuyse and Jones, 1996; Dodds et al., 
2002). 

Our model results to predict RCHL biomass were 
low (R 2 = 0.0 to 0.08), but similar to those determined 
in other agricultural streams (R 2 = 0.03) by Munn 
et al. (2010). Other studies have found stronger rela- 
tions with nutrients and RCHL biomass with R 2 
ranging from about 0.11 to 0.60 (Lohman et al., 1992; 
Chetelat et al., 1999; Dodds et al., 2006). Munn et al. 
(2010) found a negative relation between nutrients 
and RCHL in agricultural streams with open canopy, 
suggesting that periphyton can reduce TN and TP 
concentrations. However, performance of our model 
did not improve when open vs. closed canopy sites 
were analyzed separately. 

One counterintuitive result is that model perfor- 
mance generally did not improve when applied to 
individual regions. This is contrary to Stevenson 
et al. (2006) who found correlations between mea- 
sures of algal biomass and nutrients were higher 
when observations were constrained to a specific 
region vs. when data were combined. Perhaps this 
difference may have to do with our site selection that 
targeted agricultural watersheds where streams may 
have similar conditions including elevated nutrients 
and altered habitat. For example, we found no statis- 
tical differences among regions in some key environ- 
mental variables including conductivity, turbidity, 
and percent canopy all of which have been shown to 
be related to algal biomass in other streams (Biggs, 
1995; Figueroa-Nieves etal, 2006; Munn et al, 
2010). Highest mean concentrations of TN and TP 
were observed in the UMIS. However, these sites had 
the lowest mean RCHL compared with other regions, 
which may be due to generally finer substrate in 
UMIS (Biggs and Shand, 1987). 

Our best model integrated all three biotic mea- 
sures into a comprehensive index of stream eutrophi- 
cation. EI was related to nutrient concentrations (TN 
and TP) (Figure 6), but nutrient concentrations alone 
provide a poor explanation of aquatic plant growth 
responses in these streams. Adding environmental 
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factors that related to light availability (percent can- 
opy), streamflow conditions during sampling (median 
depth), and recent streamflow variability (daily 
change in stage) improve the model estimates of 
eutrophic conditions accounting for 50% of the vari- 
ance in EI scores. Aquatic plant growth increased 
with TN and decreased with percent canopy as 
expected, but it increased with median depth and 
daily change in stage. Depth appears to be a surro- 
gate for stream size, whereas daily change in stage is 
a measure of water motion that has been shown to 
enhance nutrient uptake and availability up to a 
point when the force of moving water is too great and 
plant growth is sheared from the substrate (Horner 
et al., 1990; Chambers et al., 1991). 

The findings of this study suggest that trophic- 
state classification of streams is more appropriately 
based on areal plant growth than on nutrient concen- 
trations of TN and TP. This would be especially true 
for streams that contain an abundance of macro- 
phytes. Future studies may want to consider separate 
visual coverage of macrophyte and macroalgae 
because rooted macrophytes generally uptake nutri- 
ents from bed sediments through their roots whereas 
macroalgae uptake nutrients from the water column 
(Chambers et al. , 1999). Each type of response 
(RCHL, SCHL, and AQM) may have distinct controls, 
but overall these responses are regulated by nutri- 
ents, light, and streamflow. Nutrient criteria may be 
justified as a management tool, but do not serve as a 
precise basis for classifying the eutrophic status of 
streams. 



SUPPORTING INFORMATION 

Additional Supporting Information may be found 
in the online version of this article: 

Appendix SI. Stream name, study unit, USGS site 
number, site location, sample year, and stage device 
used (G = USGS gage or S = stage recorder), n = 70. 

Please note: Neither AWRA nor Wiley-Blackwell is 
responsible for the content or functionality of any 
supporting materials supplied by the authors. Any 
queries (other than missing material) should be direc- 
ted to the corresponding author for the article. 



ACKNOWLEDGMENTS 

We thank the many hydrologists and biologists with the 
NAWQA Program who collected data used in this report. Mark 
Munn and Terry Short provided guidance in site selection and data 
interpretation. Christopher Mebane, Ross Dickinson, and Dorene 
MacCoy assisted with data collection, compilation, and graphics. 



Kathy Lee, Billy Justus, Jim Petersen, and Suzanne Femmer were 
instrumental in collecting data and provided valuable information 
about their respective study areas. Jill Frankforter provided valu- 
able data checking and summaries. Bob Black assisted with piece- 
wise regression analysis. Colleague reviews by Billy Justus, 
Christopher Mebane, and Mark Munn improved the quality of the 
report. Three anonymous reviewers and the associate editor are 
gratefully acknowledged for their suggested improvements to the 
manuscript. Funding was provided by the U.S. Geological Survey, 
National Water Quality Assessment Program. Any use of trade, 
product, or firm names in this publication is for descriptive pur- 
poses only and does not imply endorsement by the U.S. Govern- 
ment. 



LITERATURE CITED 

Allan, J.D., 1995. Stream Ecology: Structure and Function of Run- 
ning Waters. Kluwer Academic Publishers, Dordrect, The Neth- 
erlands. 

Arar, E.J. and G.B. Colling, 1997. U.S. Environmental Protection 
Agency Method 445.0, In Vitro Determination of Chlorophyll a 
and Pheophytin a in Marine and Freshwater Algae by Fluores- 
cence, Revision 1.2. U.S. Environmental Protection Agency, 
National Exposure Research Laboratory, Office of Research and 
Development, Cincinnati, Ohio, http://www.epa.gov/microbes/ 
m445_0.pdf, accessed February 2010. 

Biggs, B.J.F., 1995. The Contribution of Flood Disturbance, Catch- 
ment Geology, and Land Use to the Habitat Template of Periph- 
yton in Stream Ecosystems. Freshwater Biology 33:419-438. 

Biggs, B.J.F., 1996a. Patterns in Benthic Algae in Streams. In: 
Algal Ecology: Freshwater Benthic Ecosystems, J.R. Stevenson, 
M.L. Bothwell, and R.L. Lowe (Editors). Academic Press, San 
Diego, pp. 31-76. 

Biggs, B.J.F., 1996b. Hydraulic Habitat of Plants in Streams. Regu- 
lated Rivers: Research and Management 12:131-144. 

Biggs, B.J.F. and M.E. Close, 1989. Periphyton Biomass Dynamics 
in Gravel Bed Rivers: The Relative Effects of Flows and Nutri- 
ents. Freshwater Biology 22:209-231. 

Biggs, B.J.F. and B.I. Shand, 1987. Biological Communities and 
the Potential Effects of Power Development in the Lower Clutha 
River, Otago. Ministry of Works and Development, Christ- 
church, New Zealand. 

Borchardt, M.A., 1996. Nutrients. In: Algal Ecology: Freshwater 
Benthic Ecosystems, R.J. Stevenson, M.L. Bothwell, and R.L. 
Lowe (Editors). Academic Press, San Diego, pp. 183-227. 

Brightbill, R.A. and M.D. Munn, 2008. Environmental and Biologi- 
cal Data of the Nutrient Enrichment Effects on Stream Ecosys- 
tems Project of the National Water Quality Assessment 
Program, 2003-04. U.S. Geological Survey Data Series 345, 
Tacoma, Washington, 12 pp. http://pubs.usgs.gov/ds/345, acces- 
sed February 2010. 

Carpenter, K.D. and I.R. Waite, 2000. Relations of Habitat-Specific 
Algal Assemblages to Land Use and Water Chemistry in the 
Willamette Basin, Oregon. Environmental Monitoring and 
Assessment 64:247-257. 

Carpenter, S.R. and D.M. Lodge, 1986. Effects of Submersed 
Macrophytes on Ecosystem Processes. Aquatic Botany 26:341- 
370. 

Cattaneo, A. and M.C. Amireault, 1992. How Artificial Are Artifi- 
cial Substrata for Periphyton? Journal of the North American 
Benthological Society ll(2):244-256. 

Chambers, PA., R.E. DeWreede, EA. Irlandi, and H. Vandermeu- 
len, 1999. Management Issues in Aquatic Macrophyte Ecology: 
A Canadian Perspective. Canadian Journal of Botany 77(4):471- 
487. 



Journal of the American Water Resources Association 



511 



JAWRA 



Maret, Konrad, and Tranmer 



Chambers, P.A., E.E. Prepas, H.R. Hamilton, and M.L. Bothwell, 
1991. Current Velocity and Its Effect on Aquatic Macrophytes 
in Flowing Waters. Ecological Applications l(3):249-257. 

Chetelat, J., F.R. Pick, A. Morin, and P.B. Hamilton, 1999. Periph- 
yton Biomass and Community Composition in Rivers of Differ- 
ent Nutrient Status. Canadian Journal of Fisheries and Aquatic 
Sciences 56:560-569. 

Clark, G.M., D.K. Mueller, and M.A. Mast, 2000. Nutrient Concen- 
trations and Yields in Undeveloped Stream Basins of the United 
States. Journal of the American Water Resources Association 
36(4):849-860. 

Cuffney, T.F., 2003. User's Manual for the National Water-Quality 
Assessment Program Invertebrate Data Analysis System (IDAS) 
Software: Version 3. U.S. Geological Survey, Raleigh, North 
Carolina, 103 pp. Open-File Report 03-172, http://nc.water.usgs. 
gov/reports/ofr03172, accessed February 2010. 

Dodds, W.K., W.W. Bouska, J.L. Eitzmann, T.J. Pilger, K.L. 
Pitts, A.J. Riley, J.T. Scholesser, and D.J. Thornbrugh, 2009. 
Eutrophication of U.S. Freshwaters: Analysis of Potential Eco- 
nomic Damages. Environmental Science and Technology 
43(1):12-19. 

Dodds, W.K., J.R. Jones, and E.B. Welch, 1998. Suggested Classifi- 
cation of Stream Trophic State: Distributions of Temperate 
Stream Types by Chlorophyll, Total Nitrogen, and Phosphorous. 
Water Research 32:1455-1462. 

Dodds, W.K., V.H. Smith, and K. Lohman, 2002. Nitrogen and 
Phosphorus Relationships to Benthic Algal Biomass in Temper- 
ate Streams. Canadian Journal of Fisheries and Aquatic Sci- 
ences 59:865-874. 

Dodds, W.K., V.H. Smith, and K. Lohman, 2006. Erratum: Nitro- 
gen and Phosphorus Relationships to Benthic Algal Biomass in 
Temperate Streams. Canadian Journal of Fisheries and Aquatic 
Sciences 63:1190-1191. 

Dodds, W.K. and E.B. Welch, 2000. Establishing Nutrient Criteria 
in Streams. Journal of the North American Benthological Soci- 
ety 19(1):186-196. 

Figueroa-Nieves, D., TV. Royer, and M.B. David, 2006. Controls on 
Chlorophyll-a in Nutrient-Rich Agricultural Streams in Illinois, 
USA. Hydrobiologia 568:287-298. 

Fitzpatrick, FA., I.R. Waite, P.J. DArconte, M.R. Meador, M.A. 
Maupin, and M.E. Gurtz, 1998. Revised Methods for Character- 
izing Stream Habitat in the National Water-Quality Assessment 
Program. Water-Resource Investigation, Reston, Virginia, 67 pp. 
Report 98-4052, http://pubs.usgs.gov/wri/wri984052/index.html, 
accessed February 2010. 

Francoeur, S.N., B.J.F. Biggs, R.A. Smith, and R.L. Lowe, 1999. 
Nutrient Limitation of Algal Biomass Accrual in Streams: Sea- 
sonal Patterns and a Comparison of Methods. Journal of the 
North American Benthological Society 18(2):242-260. 

Helsel, D.R. and R.M. Hirsch, 2002. Statistical Methods in Water 
Resources, Techniques of Water-Resources Investigations of the 
United States Geological Survey, Book 4, Hydrologic Analysis 
and Interpretation. http://water.usgs.gov/pubs/twri/twri4a3/, 
accessed February 2010. 

Hering, D., R.K. Johnson, S. Kramm, S. Schmutz, K. Szoszkiewicz, 
and P.F.M. Verdonschot, 2006. Assessment of European 
Streams With Diatoms, Macrophytes, Macroinvertebrates and 
Fish: A Comparative Metric-Based Analysis of Organism Re- 
ponse to Stress. Freshwater Biology 51:1757-1785. 

Horner, R.R., E.B. Welch, M.R. Seeley, and J.M. Jacoby, 1990. 
Responses of Periphyton to Changes in Current Velocity, Sus- 
pended Sediment and Phosphorus Concentration. Freshwater 
Biology 24(2):215-232. 

Hughes, R.M., P.R. Kaufmann, A.T. Helihy, T.M. Kincaid, L. Rey- 
nolds, and D.P. Larsen, 1998. A Process for Developing and 
Evaluating Indices of Fish Assemblage Integrity. Canadian 
Journal of Fisheries and Aquatic Sciences 55:1618-1631. 



Jacoby, J.M., 1987. Alterations in Periphyton Characteristics Due 
to Grazing in a Cascade Foothill Stream. Freshwater Biology 
18:495-508. 

Koetsier, P., 2005. The Effects of Disturbance Time Interval on 
Algal Biomass in a Small Idaho Stream. Northwest Science 
79(4):211-216. 

Konrad, CP., D.B. Booth, S.J. Burges, and D.R. Montgomery, 
2002. Partial Entrainment of Gravel Bars During Floods. Water 
Resources Research 38(7):1104. 

Lamberti, G.A. and V.H. Resh, 1983. Stream Periphyton and Insect 
Herbivores: An Experimental Study of Grazing by a Caddisfly 
Population. Ecology 64(5):1124-1135. 

Lohman, K., J.R. Jones, and B.D. Perkins, 1992. Effects of Nutrient 
Enrichment and Flood Frequency on Periphyton Biomass in 
Northern Ozark Streams. Canadian Journal of Fisheries and 
Aquatic Sciences 49(6):1198-1205. 

Maret, T.R., D.E. MacCoy, and D.M. Carlisle, 2008. Long-Term 
Water Quality and Biological Responses to Multiple Best Man- 
agement Practices in Rock Creek, Idaho. Journal of the Ameri- 
can Water Resources Association 44(5):1248-1269. 

McMahon, G., J.D. Bales, J.F. Coles, E.M.P. Giddings, and H. Zap- 
pia, 2003. Use of Stage Data to Characterize Hydrologic Condi- 
tions in an Urbanizing Environment. Journal of the American 
Water Resources Association 39(6):1529-1546. 

Morin, A. and A. Cattaneo, 1992. Factors Affecting Sampling Vari- 
ability of Freshwater Periphyton and the Power of Periphyton 
Studies. Canadian Journal of Fisheries and Aquatic Sciences 
49(8):1695-1703. 

Moulton, S.R., II, J.L. Carter, S.A. Grotheer, T.F. Cuffney, and 
T.M. Short, 2000. Methods for Analysis by the U.S. Geological 
Survey National Water-Quality Assessment Laboratory - 
Processing, Taxonomy, and Quality Control of Benthic Macro- 
invertebrate Samples. U.S. Geological Survey, Raleigh, North 
Carolina, 49 pp. Open-File Report 00-212, http://nwql.usgs.gov/ 
Public/pubs/OFR00-212.html, accessed February 2010. 

Moulton, S.R., II, J.G. Kennen, R.M. Goldstein, and JA. Ham- 
brook, 2002. Revised Protocols for Sampling Algal, Inverte- 
brate, and Fish Communities as Part of the National Water- 
Quality Assessment Program. U.S. Geological Survey, Raleigh, 
North Carolina, 75 pp. Open-File Report 02-150, http:// 
pubs.usgs.gov/of/2002/ofr-02-150/index.html, accessed February 
2010. 

Mueller, D.K. and N.E. Spahr, 2006. Nutrients in Streams and 
Rivers Across the Nation - 1992-2001. U.S. Geological Survey, 
Denver, Colorado, 44 pp. Scientific Investigations Report 2006- 
5107, http://pubs.er.usgs.gov/usgspubs/sir/sir20065107, accessed 
February 2010. 

Munn, M.D., R.W. Black, and S.J. Gruber, 2002. Response of Ben- 
thic Algae to Environmental Gradients in an Agriculturally 
Dominated Landscape. Journal of the North American Bentho- 
logical Society 21(2):221-237. 

Munn, M.D., J.W. Frey, and A.J. Tesoriero, 2010. The Influence of 
Nutrients and Physical Habitat in Regulating Algal Biomass in 
Agricultural Streams. Environmental Management, doi: 
10. 1007/s00267-010-9435-0. 

Munn, M.D., L.L. Osborne, and M.J. Wiley, 1989. Factors Influenc- 
ing Periphyton Growth in Agricultural Streams of Central Illi- 
nois. Hydrobiologia 174:89-97. 

Munn, M.D., I.R. Waite, D.P. Larsen, and AT. Herlihy, 2009. The 
Relative Influence of Geographic Location and Reach-Scale Hab- 
itat on Benthic Invertebrate Assemblages in Six Ecoregions. 
Environmental Monitoring and Assessment 154:1-14. 

Myers, R., J. Harrison, S.K. Parkinson, B. Hoelscher, J. Naymik, 
and S.E. Parkinson, 2003. Pollutant Transport and Processing 
in the Hells Canyon Complex. Idaho Power Company, Boise, 
Idaho, 218 pp., Technical Report Appendix E.2.2-2, Hells Can- 
yon Hydroelectric Project New License Application (FERC No. 



JAWRA 



512 



Journal of the American Water Resources Association 



Influence of Environmental Factors on Biotic Responses to Nutrient Enrichment in Agricultural Streams 



1971). http://www.idahopower.com/pdfs/Relicensing/hellscanyon/ 
hellspdfs/TechAppendices/Water%20Quality/e22_02.pdf, accessed 
February 2010. 

van Nieuwenhuyse, E.E. and J.R. Jones, 1996. Phosphorus-Chloro- 
phyll Relationship in Temperate Streams and Its Variation 
With Stream Catchment Area. Canadian Journal of Fisheries 
and Aquatic Sciences 53:99-105. 

Omernik, J.M., 1987. Map Supplement: Ecoregions of the Conter- 
minous United States. Scale 1:750,000. Annals of the Associa- 
tion of American Geographers 77:118-125. 

Patton, C.J. and J.R. Kryskalla, 2003. Methods of Analysis by the 
U.S. Geological Survey National Water Quality Laboratory - 
Evaluation of Alkaline Persulfate Digestion as an Alternative to 
Kjeldahl Digestion for Determination of Total and Dissolved 
Nitrogen and Phosphorus in Water. U.S. Geological Survey, 
Open-File Report 03-4174. Denver, Colorado, 33 pp. http:// 
nwql.usgs.gov/pubs-wrir.shtml, accessed February 2010. 

Patton, C.J. and E.P. Truitt, 2000. Methods of Analysis by the U.S. 
Geological Survey National Water Quality Laboratory-Determi- 
nation of Ammonium Plus Organic Nitrogen by a Kjeldahl 
Digestion Method and an Automated Photometric Finish That 
Includes Digest Cleanup by Gas Diffusion. U.S. Geological Sur- 
vey, Open-File Report 00-170. Denver, Colorado, 31 pp. http:// 
nwql.usgs.gov/WRIR-03-4174.shtml, accessed February 2010. 

Perkins, B.D. and J.R. Jones, 1994. Temporal Variability in a Mid- 
western Stream During Spring. Verhandlungen des Interna- 
tional en Verein Limnologie 25:1471-1476. 

Poff, N.L., 1997. Landscape Filters and Species Traits: Towards 
Mechanistic Understanding and Prediction in Stream Ecology. 
Journal of the North American Benthological Society 16(2):391- 
409. 

Porter, S.D., D.K. Mueller, N.E. Spahr, M.D. Munn, and N.M. Dub- 
rovsky, 2008. Efficacy of Algal Metrics for Assessing Nutrient 
and Organic Enrichment in Flowing Waters. Freshwater Biol- 
ogy 53(5):1036-1054. 

R Development Core Team, 2005. R: A Language and Environment 
for Statistical Computing. R Foundation for Statistical Comput- 
ing, Vienna. ISBN 3-900051-07-0. http://www.r-project.org/, 
accessed February 2010. 

Rantz, S.E., et al., 1982. Measurement and Computation of Stream- 
flow: Volume 1. Measurement of Stage and Discharge. U.S. Geo- 
logical Survey Water-Supply Paper 2175, Reston, Virginia, 313 
pp. http://pubs.usgs.gov/wsp/wsp2175/, accessed February 2010. 

Shelton, L.R., 1994. Field Guide for Collecting and Processing 
Stream-Water Samples for the National Water-Quality Assess- 
ment Program. U.S. Geological Survey Open-File Report 94-455. 
U.S. Geological Survey, Reston, Virginia, 42 pp. http://pubs.er. 
usgs.gov/usgspubs/ofr/ofr94455, accessed February 2010. 

Snyder, E.B., C.T. Robinson, GW. Minshall, and S.R. Rushforth, 
2002. Regional Patterns in Periphyton Accrual and Diatom 
Assemblage Structure in a Heterogeneous Nutrient Landscape. 
Canadian Journal of Fisheries and Aquatic Sciences 59:564-577. 

Stanley, E.H., S.G. Fisher, and N.B. Grimm, 1997. Ecosystem 
Expansion and Contraction in Streams. Bioscience 47:427-435. 

Stevenson, R.J., S.T. Rier, CM. Riseng, R.E. Schultz, and M.J. Wi- 
ley, 2006. Comparing Effects of Nutrients on Algal Biomass in 
Streams in Two Regions With Different Disturbance Regimes 
and With Applications for Developing Nutrient Criteria. Hydro- 
biologia 561:149-165. 

Sumner, W.T. and CD. Mclntire, 1982. Grazer-Periphyton Interac- 
tions in Laboratory Streams. Archiv fur Hydrobiologie 93:135- 
157. 

Suplee, M.W., V. Watson, M. Teply, and H. McKee, 2009. How 
Green Is Too Green? Public Opinion of What Constitutes 
Undesirable Algae Levels in Streams. Journal of the American 
Water Resources Association 45(1):123-140. 



Toms, J.D. and M.L. Lesperance, 2003. Piecewise Regression: A 
Tool for Identifying Ecological Thresholds. Ecology 84:2034- 
2041. 

USEPA (U.S. Environmental Protection Agency), 2000a. Nutrient 
Criteria Technical Guidance Manual, River and Streams. Office 
of Water, EPA-822-B-00-002, Washington, D.C http://www.epa. 
gov/waterscience/criteria/nutrient/guidance/rivers, accessed Feb- 
ruary 2010. 

USEPA (U.S. Environmental Protection Agency), 2000b. Ecoregion- 
al Nutrient Criteria for Rivers and Streams. Office of Water, 
EPA 822-B-00-016, EPA 822-B-00-018, EPA 822-B-00-020, 
Washington, D.C. http://www.epa.gov/waterscience/criteria/ 
nutrient/ecoregions/rivers/index.html, accessed February 2010. 

USEPA (U.S. Environmental Protection Agency), 2006. Wadeable 
Streams Assessment. Office of Water, EPA 841-B-06-002, Wash- 
ington, D.C. http://www.epa.gov/owow/streamsurvey/, accessed 
February 2010. 

Wang, L., D.M. Robertson, and P.J. Garrison, 2007. Linkages 
Between Nutrients and Assemblages of Macroinvertebrates and 
Fish in Wadeable Streams: Implication to Nutrient Criteria 
Development. Environmental Management 39(2):194-212. 

Watzin, M.C and A.W. Mcintosh, 1999. Aquatic Ecosystems in 
Agricultural Landscapes: A Review of Ecological Indicators and 
Achievable Ecological Outcomes. Journal of Soil and Water Con- 
servation 54(4):636-644. 

Welch, E.B., 1980. Ecological Effects of Wastewater. Cambridge 
University Press, New York, 337 pp. 

Wilkinson, L., 2004. SYSTAT for Windows - Statistics, Version 
11.0. SYSTAT, Evanston, Illinois. 

Wolman, M.G., 1954. A Method for Sampling Coarse River-Bed 
Material. Transactions of the American Geophysical Union 
35:951-956. 

Wolock, D.M., 2003. Base-Flow Index Grid for the Conterminous 
United States. U.S. Geological Survey, Reston, Virginia, http:// 
water.usgs.gov/GIS/metadata/usgswrd/XML/bfi48grd.xml, access- 
ed February 2010. 



Journal of the American Water Resources Association 



513 



JAWRA 



